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(N 
O 

o 

(N 



> 
in 

m 

o 

CN 
O 

S3 

Of 

6 

CO 

C3 



X 



R. Napiwotzki 1 , D. Koester 2 , G. Nelemans 3 , L. Yungelson 4 , N. Christlieb 5 , A. Renzini 6 , D. Reimers 5 

H. Drechsel 1 , and B. Leibundgut 6 



Dr. Remeis-Sternwarte, Astronomisches Institut der Universitat Erlangen-Niirnberg, Sternwartstr. 7, D-96049 
Bamberg, Germany; e-mail: napiwotzki@sternwarte.uni-erlangen.de 

Institut fur Theoretische Physik und Astrophysik, Universitat Kiel, D-24098 Kiel, Germany 
Institute of Astronomy, Madingley Road, CB3 OHA, Cambridge, UK 

Institute of Astronomy of the Russian Academy of Sciences, 48 Pyatnitskaya Str., 109017 Moscow, Russia 
Hamburger Sternwarte, Universitat Hamburg, Gojenbergsweg 112, D-21029 Hamburg, Germany 
European Southern Observatory, Karl-Schwarzschild-Str. 2, D-85748 Garching, Germany 



Abstract. In the course of our search for double degenerate (DD) binaries as potential progenitors of type la 
supernovae with the UVES spectrograph at the ESO VLT (ESO SN la Progenitor surveY - SPY) we discovered 
that the white dwarf HE 1414-0848 is a double-lined DA+DA binary with an orbital period of P = 12 h 25 m 39 3 . 
Semi-amplitudes of 128 kms -1 and 100 kms -1 are derived for the individual components. The amplitude ratio 
and the measured difference in gravitational redshift is used to estimate the masses of the individual compo- 
nents: O.55M and O.71M . Hence the total mass of the HE 1414-0848 system is 1.26M , only 10% below the 
Chandrasekhar limit. The results of a model atmosphere analysis are consistent with our mass estimated from the 
orbit. Temperatures of the individual components are also determined. Possible scenarios for the formation of this 
system are discussed. The system will merge due to loss of angular momentum via gravitational wave radiation 
after two Hubble times. HE 1414— 0848 does not qualify as a SNIa progenitor, but it is the most massive close 
DD known today. 
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1. Introduction 

Supernovae of type la (SNIa) play an outstanding role 
in the study of cosmic evolution. In particular they are 
regarded as one of the best standard candles to determine 
the cosmological parameters H , fl and A (e.g. Riess et 
al. 199S, Leibundgut 2001). However, the nature of their 
progenitors remains a mystery (e.g. Livio 200C). 

There is general consensus that the SN event is due 
to the thermonuclear explosion of a white dwarf (WD) 
when a critical mass (likely the Chandrasekhar limit) is 
reached, but the nature of the progenitor system remains 
unclear. While it must be a binary, with matter being 
transferred to the WD from a companion until the critical 
mass is reached, two main options exist for the nature 
of the mass donator: either another WD in the so-called 



or a red giant/subgiant in the so-cal led single degenerate 
(SD) scenario (Whelan & Iben |l973|) . 



double degenerate (DD) scenario (Iben & Tutukov |1984| ) 



Send offprint requests to: R. Napiwotzki 

* Based on data obtained at the Paranal Observatory of the 
European Southern Observatory for programs No. 165.H-0588, 
266.D-5658, and 167.D-0407 



As progenitor candidate, the DD model considers a bi- 
nary with the total mass of the white dwarf components 
larger than the Chandrasekhar mass, which merges in less 
than a Hubble time due to the loss of angular momen- 
tum via gravitational wave radiation. Several systematic 
radial velocity (RV) searches for DDs have already been 
undertaken starting in the mid 1980's (see Marsh 2000 
for a review). By now, combining all the surveys « 180 
WDs have been checked for RV variations yielding a sam- 
ple of 18 DDs with P < 6 d 3 (Marsh |2000| Maxted et al. 
|2000a ). Six of the 18 known DDs are double-lined sys- 
tems (in three of these systems the companion is barely 
detectable). None of the 18 known DD systems seems mas- 
sive enough to qualify as a SN la precursor. This is not 
surprising, as theoretical simulations suggest that only a 
few percent of all close DDs are potenti al SN Ia progen- 
itors (Iben et al. 



1997 



Nelemans et al. p00l|) . Recently, 
the binary KPD 1930+2752, consisting of a subluminous 
B (sdB) star and an invisible white dwarf component, 
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was proposed as potential SN la progenitor by Maxted et 
al. (20001). The system mass exceeds the Chandrasekhar 
limit, and the system will merge within a Hubble time. 
However, it is not clear, if the merger event will produce 
a SNIa explosion (Ergma et al. 2001). 

In order to perform a definite test of the DD scenario 
we have embarked on a large spectroscopic survey of 1500 
white dwarfs using the UVES spectrograph at the ESO 
VLT UT2 (Kueyen) to search for white dwarfs and pre- 
white dwarfs with variable RVs (ESO SNIa Progenitor 
surveY - SPY). The SPY project yields a wealth of new 
RV variable DDs (Napiwotzki et al. 2001b| ). An analysis 
of the subdwarf B + white dwarf system HE 1047-0436 
was p resente d in the first paper of this series (Napiwotzki 
et al. bOOlal). 



Among the newly detected DDs are six new double- 
lined binaries. In all systems both WD components can 
be easily recognised. Here we report on follow-up spec- 
troscopy of HE 1414-0848 (cf. Sect. 2). The determina- 
tion of the radial velocity curves, orbital parameters, and 
masses for both components is described in Sect. 3. A 
model atmosphere analysis is carried out and the results 
are compared with the analysis of the RV curves in Sect. 4. 
We finish with a discussion of the HE 1414—0848 system 
and possible formation scenarios. 



2. Observations and Data Analysis 

HE 1414-0848 (a 20 oo = 14 h 16 m 52.1 s , <S 2000 = -9°02'04", 
B pg =16™2) was disc overed by the Hambur g ES O survey 
(HES; Wisotzki et al. |2000| Christlieb et al. $00l\ ) as a po- 
tential cool white dwarf and, therefore, was included in our 
survey. The survey spectra showed that it is a double- lined 
system consisting of two DA white dwarfs. Both compo- 
nents were visible in the H Q core separated by 4.6 A, cor- 
responding to 210km/s, in both discovery spectra. This 
indicated high orbital velocity and a short period, making 
HE 1414—0848 a high-priority target for follow-up studies. 

Fifteen high resolution (0.3A) echelle spectra of 
HE 1414-0848 have been secured with VLT-UVES be- 
tween March 9 and April 5, 2001 in service mode. Details 
on the observational set up of the UVES instrument and 
the data reduction can be found in Koester et al. (2001) 
and Napiwotzki et al. ( 2001b[ ). The chosen observing strat- 
egy was to start with several closely spaced exposures dur- 
ing the first night. This should provide us with a rough 
first estimate of the period and radial velocity curve, which 
could be used to extrapolate the orbit to the next nights. 
Consequently, a decreasing number of spectra was taken 
during the following nights, finishing with single exposures 
separated by more than one day. This strategy in conjunc- 
tion with the high quality of the UVES data allowed the 
determination of accurate orbital parameters with a rela- 
tively small number of spectra. 

Six H Q spectra taken during the first night of the 
follow-up observations are displayed in Fig. |lj together 
with a fit of both components as described in Sect. 0. 
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(*= 0.459) 
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1977.810 
(*= 0.737) 



1977.876 
($= 0.864) 
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Fig. 1. H Q region of the spectra covering a timespan of 
5 h in the night of March 8/9, 2001 together with the fits 
described in Sect. ||. The numbers indicate the Julian date 
of the exposures and the orbital phase (f> (cf. Fig. ||). Star 
"A" in the text corresponds to the blueward shifted line 
core in, e.g., the JD 2451977.773 exposure, "B" to the red 
component. Observed spectra were smoothed with a three 
pixel boxcar (0.09 A). 



The rapid change of the spectral appearance due to the 
orbital motion is obvious. 



3. Radial velocity curve and orbital parameters 

Although the line cores of both components are similar, 
the left one in Fig. [I] is slightly deeper and broader (called 
component A further on). We used this in a first step to 
identify both stars in each spectrum. In a second step we 
fitted the central region (±20 A) by using two Gaussians 
(one for every component) , a Lorentzian to model the line 
wings (parameters of the Lorentzian were determined from 
a fit covering ±50 A and then held fixed), and a linear 
polynomial to reproduce the overall spectral trend. Fitting 
of the H Q profiles was performed with a downhill simplex 
algorithm, which will be described in a forthcoming paper. 

In a third step we fitted RV curves for a range of pe- 
riods and produced a "power spectrum" indicating the 
fit quality as a function of period (like the one presented 
in Fig. ||). In a fourth step we inspected the phased RV 
curves using periods corresponding to peaks in the power 
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Table 1. Fitted heliocentric radial velocities. Hel.JD was 
calculated for the middle of the exposures. 



Hel.JD 


Heliocentric RV 


-2,400,000 


A 


B 


51684.68664 


-110±4 


98±4 


51687.61061 


129±4 


-80±4 


51977.66606 


47±9 


-3±11 


51977.67313 


24±4 


-4±6 


51977.73344 


-59±5 


70±5 


51977.77310 


-106±4 


92±5 


51977.81012 


-119±9 


104±11 


51977.87602 


-98±5 


91±6 


51978.67757 


67±4 


-34±6 


51978.75510 


-50±4 


56±4 


51978.85895 


-124±5 


112±5 


51979.68326 


92±8 


-64±9 


51979.88119 


-123±7 


113±8 


51987.70636 


-95±3 


92±5 


52003.59657 


-29±3 


39±4 


52003.82394 


-25±3 


44±3 


52004.58719 


41±4 


-20±5 


52137.53531 


127±3 


-87±4 




0.5 0-5 1 1.5 

Phase 

Fig. 3. Measured radial velocities as a function of orbital 
phase and fitted sine curves for HE 1414-0848. Filled cir- 
cles correspond to component "A" and open rectangles 
to "B". Note the difference of the "systemic velocities" 
7o between both components caused by the gravitational 
redshift. 



Table 2. Orbital parameters of the HE 1414-0848 sys- 
tem. Given are the "system" velocities 70, RV amplitudes 
Ka/b, the x 2 values of the fits and the number of data 
points N. 



Comp. 


70 [km/s] 


Ka/b [km/s] 


N 


x 2 


A 


2.1 ± 1.1 


128.2 ± 1.6 


18 


23.7 


B 


16.4 ± 1.3 


99.9 ±1.7 


18 


25.2 





period/ days 



Fig. 2. Power spectrum of the RV measurements of com- 
ponent A. The inset shows details of the region around 
the main peak. 



spectrum and compared them with the measured values. 
Most could be ruled out, because of discrepant RV mea- 
surements. Finally, one unambiguous solution remained. 

In the fifth and last step we extrapolated the de- 
rived RV curves back to the time of the discovery ob- 
servations and to an additional spectrum taken in August 
2001. We used the phase information to identify the com- 



ponents in these exposures, and added the RVs to our 
analysis. Individual measurements are listed in Table 
Since the spectra span more than one year, a very accu- 
rate period could be determined: P = 12 h 25 m 39 s . Results 
for both components agree within 1 s. Semi-amplitudes 
Ka/b and the "system velocities" 70 are given in Table |[ 
Accordingly the ephemcris for the time To defined as the 
conjunction time at which star A moves from the blue side 
to the red side of the RV curve (i.e. star A is closest to 
the observer) is 



Hel.JD (T ) = 2451976.911 + 0.51781 xE 
±0.001 ±0.00001 



(1) 



Error limits are given in the second line. The mass ratio 
of both white dwarfs can be computed from the ratio of 
RV amplitude: 



= 1.28 ±0.02. 



(2) 



M B = Ka 
M a ~ K B 

From Fig. |3] and Table || it is evident that the "system 
velocities" 70 derived for component A and B differ much 
more than naively expected considering the error bars: 
70(B) — 70(A) = 14.3 ± 1.7km/s . However, this is easily 
explained by the mass dependent gravitational redshift of 
white dwarfs 

GM 



Rc 2 



(3) 



This offers the opportunity to determine masses of the 
individual white dwarfs in the HE 1414—0848 binary. For 
a given mass-radius relation gravitational redshift can be 
computed as a function of white dwarf mass. Since the 
mass ratio is given by Eq. [2J only one combination of 
masses can fulfil both constraints. Since HE 1414—0848 
is the outcome of common envelope evolution, it is not 
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clear which mass-radius relation to use. We performed 
the analysis with mass-radius relations constructed for 
T cff = 9800 K (cf. Sect, g) from three different white 
dwarf cooling tracks, which should give us an estimate 
of systematic errors introduced by a particular choice. 
We used the tracks computed by Blocker et al. (1997), 
which result from the self-consistent calculation of sin- 
gle star evolution and cooling curves computed by Wood 
(1995) for white dwarfs with a "thick" hydrogen layer 
(M H = 10" 4 M WD ) and a "thin" layer (M H = 0). The re- 
sulting white dwarf masses are Ma = 0.55 ± 0.03M Q and 
M B = 0.71 ± 0.03M Q for Wood's models with (M H = 0). 
The resulting masses are 0.01M Q larger/lower if we adopt 
the ( Mh = 10 _4 A/wd) models or the tracks of Blocker 
et al. ( 1997 ), respectively. Thus the choice of a particu- 
lar mass-radius relation is of little importance. The mass 
sum of both white dwarfs is M — 1.26 ± 0.06M Q (er- 
rors propagated from the orbital parameters + possible 
systematic errors). Thus HE 1414— 0848 is the most mas- 
sive DD ever found with a total mass only 10% below the 
Chandrasekhar limit! 

The separation between both white dwarfs is quite 
small, only 2.9R®. Thus HE 1414-0848 obviously under- 
went phases of strong binary interaction in its history. 
From the size of the orbits and the period orbital veloc- 
ities can be computed: 164 km s _1 for component A and 
123 km s _1 for B. The comparison with the observed RV 
amplitudes allows us to determine the inclination of this 
system as i = 52°. 



[IE1414-0K4N PI.- 




AX [A] 



Fig. 4. Fit of four spectra of HE 1414-0848 taken close 
to conjunction (cf. Table |J). 



4. Spectroscopic analysis 

4.1. Spectral fitting with one model 

Further insight into the HE 1414—0848 system can be 
gained from a model atmosphere analysis. Since this sys- 
tem is double-lined the spectra are a superposition of both 
individual white dwarf spectra. A deconvolution is beyond 
the scope of this paper, but an analysis of the combined 
spectrum will already allow us to derive constraints on this 
system. We used four spectra obtained close to conjunc- 
tion for this purpose (Table ||) . Radial velocity differences 
are small and do not cause significant artificial broadening 
of the line profiles. 

The spectra were rebinned to a resolution of approx- 
imately 1 A. A large grid of LTE DA models computed 
with a model atmosphere code described in Finley et al. 
(11997D was used for the analysis. A simultaneous fit of the 
Balmer lines Ha to Hs using a \ 2 minimisation technique 



was performed. For details refer to Koester et al. ( 2001 
and references therein). The fits are shown in Fig. [I and 
individual results are given in Table ^. The average param- 
eters are T cS = 9800 ± 90 K and \ogg = 8. 17 ±0.09 (errors 
correspond to the scatter of individual measurements). 

Since the observed spectra are a combination of spec- 
tra of both components, the interpretation of this result 
is not straightforward. However, we expect that the anal- 
ysis results represent some sort of average of the individ- 



Table 3. Orbital phases 
displayed in Fig. [|. 



and fit results for the spectra 



spectrum 


4> 


r e a [K] 


log g [cm/s 2 ] 


PI 


0.459 


9739 ± 29 


8.12 ±0.04 


Pl_2 


0.472 


9722 ± 27 


8.21 ±0.04 


P13 


0.975 


9924 ± 26 


8.29 ±0.04 


P14 


0.449 


9822 ± 19 


8.08 ±0.02 



ual values of both components (especially, if both compo- 
nents are similar as in the case of HE 1414—0848; cf. next 



Section). We used Wood's ( 1995 ) mass radius-relation for 
(Mh = 0) to estimate a mass of 0.68 ± O.O6M0 for these 
parameters. This is consistent with the results of our anal- 
ysis of the RV curve. 

4.2. Spectral fitting with two models 

Fitting the combined spectra of the two white dwarfs can 
only give an approximate average of the atmospheric pa- 
rameters. We have therefore tried to use 14 individual 
spectra covering different phases of the orbit and fit them 
with a combination of two DA spectra (another spectrum 
- 12 - has a very perturbed Hcv profile and could not be 
used). The quality of these spectra varies and the total 
number of fit parameters is high - it is not possible to de- 
termine them unambiguously. We have therefore used the 
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Fig. 5. Combination of the best fitting model spectra for 
components A and B (T cff (A) = 8900 K and T cff (B) = 
10800 K and four different phases of the binary orbit. 

radial velocities of the individual components and their 
masses from the radial velocity curve and held them fixed 
during the fitting. Surface gravities are 8.18 and 7.87; this 
also determines the relative weight of the two model spec- 
tra from the radius, obtained from the mass-radius rela- 
tion of Wood (1995); the weight ratio is 40:60. The only 
remaining parameters are then the two effective tempera- 
tures of the models, which are obtained through a \ 2 fit - 
ting routine similar to the methods described in Koester 
et al. (2001). 

Table | gives the results for two different attempts. 
With Method 1 we have used only the inner parts of Ha 
and H/3 (80 - 100 A) to determine the atmospheric param- 
eters. In Method 2 we performed a simultaneous fit of all 
Balmer lines from Ha to H§. 

This exercise demonstrates that the individual tem- 
peratures are similar for both components and are in fact 
fairly close to the values obtained from fitting the com- 
bined spectra with one average model. However, the scat- 
ter from the different spectra is quite large - even the 
assignment, which of the two stars is the hotter one is 




log P(hr) 

Fig. 6. Period-mass distribution for the current popula- 
tion of DDs in the Galaxy for systems in which at least one 
white dwarf has a mass above 0.7M Q from the best model 



(A3) of Nelemans et al. ( 2001 ). Grey shades indicate the 
expected number of systems brighter than V = 15 m . The 
symbols show the possible mass ratios' for HE 1414—0848 
depending on which of the two components was formed 
first. 

not the same in all phases. Obviously the demand on the 
signal-to-noise ratio and on the calibration is higher for 
this method than for fitting single stars, especially if the 
two spectra are very similar as in this case. 

The most meaningful results can be expected dur- 
ing the quadrature phases ((f) — 0.25 and 0.75), when 
the RV separation between both components is largest. 
Thus not unexpectedly the scatter of the temperature de- 
terminations for the corresponding spectra P3, P4, P5, 
P6, and P10 (marked boldface in Table |^) is relatively 
small. Our best temperature estimates are 8900 ± 550 K 
for component A and 10790 ± 550 K for B from Method 2 
(9330 ± 640 K and 10310 ± 480 K if Method 1 is adopted). 
The error margins correspond to the scatter of individual 
measurements. Sample spectra fitted with these parame- 
ters for several orbital phases are displayed in Fig. |j| 

5. Discussion 



From the above derived temperatures we can estimate the 
ages of both components using cooling tracks for white 
dwarfs. Ages were interpolated from the tracks of Wood 
(1995) for the masses of the white dwarfs (temperatures 
from Method 2 adopted). Resulting ages are (1.2 ± 0.2) • 
10 9 yrs for A and (1.0±0.3) • 10 9 yrs for B. Thus the cooling 
ages of both white dwarfs are very similar and although 
B is «2000K hotter than A it could, within our error 
margins, have been formed first. 

We can compare the properties of HE 1414—0848 with 
the outcome of computations for the synthesis of th e pop - 
ulation of double white dwarfs of (Nelemans et al. 2001 ). 
In Fig. § we compare the predicted distribution over pe- 
riod and mass ratios' for systems with masses similar to 
the masses of HE 1414—0848. The mass ratio in this plot 
m/M is defined as the mass of the last formed white dwarf 
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Table 4. Results for the effective temperatures and their ler formal errors of the two DA components from spectra 
of different phases. Index (B) indicates the more massive star. Columns 3-6 correspond to Method 1, where only the 
inner parts of Ha and H/3 were used in the fit. Columns 7-10 correspond to Method 2, which used Ha to Hs. 



spectrum 


4> 


T cff (B) 


<r(B) 


T off (A) 


a(A) 


T cff (B) 


a(B) 


T cff (A) 


a(A) 


PI 


0.458 


8614 


76 


10143 


20 


10065 


78 


9427 


60 


Pl_2 


0.471 


9668 


210 


9812 


133 


9819 


204 


9735 


131 


P2 


0.581 


9950 


102 


9345 


88 


7837 


105 


9974 


12 


P3 


0.664 


9811 


148 


9756 


99 


11076 


18 


8612 


28 


P4 


0.736 


10867 


22 


8482 


41 


10916 


22 


8329 


38 


P5 


0.863 


10769 


32 


8793 


38 


9814 


113 


9798 


72 


P6 


0.847 


10006 


141 


9900 


87 


11056 


31 


9154 


24 


P7 


0.534 


9748 


171 


9858 


109 


9807 


171 


9803 


110 


P8 


0.411 


9303 


111 


10161 


48 


10198 


69 


9411 


59 


P9 


0.561 


8087 


72 


10219 


9 


10862 


20 


8574 


32 


P10 


0.761 


10100 


139 


9696 


99 


11075 


20 


8611 


31 


Pll 


0.353 


9765 


267 


9728 


172 


11087 


37 


7648 


81 


P13 


0.973 


9722 


131 


9823 


83 


10151 


62 


9683 


45 


P14 


0.447 


10005 


219 


9900 


137 


9441 


111 


10210 


48 



over the mass of the first formed white dwarf. Since in the 
case of HE 1414—0848 it is not reliable known which white 
dwarf formed first, we plotted both possible mass ratios, 
1.28 and 0.78, in Fig. |(| The population synthesis indicates 
a larger probability for the systems with m/M = 0.78. In 
Fig. ^ we give an example of a scenario leading to such a 
DD system. 

A binary with an initial semi-major axis of ~150i?© 
has components with initial masses of «4 and 3M Q . The 
primary fills its Roche lobe in the AGB stage and un- 
stable mass transfer ensues. We estimate the outcome of 
the mass transfer using an angular momentum balance 
(Nelemans et al. 200C). The orbit of the system shrinks 
slightly and the primary becomes a hot carbon-oxygen 
(CO) white dwarf. After «170 Myr the secondary fills its 
Roche lobe as an AGB star. The mass transfer is again 
unstable, the common envelope is lost, and a close binary 
consisting of two CO white dwarfs, just like HE1414-0848, 
is the result (see Fig Q). The period of the system decreases 
during the subsequent evolution due to angular momen- 
tum loss via radiation of gravitational waves. 

6. Conclusions 

We report the discovery of the double-lined radial veloc- 
ity variable binary HE 1414—0848. We measured accurate 
radial velocities for both components from high resolu- 
tion spectra and derived an orbital period of 12 h 25 m 39 s 
and semi-amplitudes of 128km/s and lOOkm/s, respec- 
tively. We combined the mass ratio computed from this 
amplitudes and the measured difference of gravitational 
redshifts and determined masses of 0.55 ± 0.03M Q and 
0.71 ± 0.03AT Q for the components using a mass-radius 
relation. This is the first application of this method to 
a DD system. The total mass of the HE 1414-0848 sys- 
tem is 1.26±0.06Mq, only 10% below the Chandrasekhar 
limit, making it the most massive DD ever found in a RV 
survey. It will merge within two Hubble times (24Gyrs). 



M,(M Q ) 



M 2 (M S ) P b (d) 



t (Myr) 




Fig. 7. A possible scenario for the formation of a system 
similar to HE1414-0848, starting with a 3.86 + 2.97M 
star in a 73.5 day orbit. The black line shows the Roche 
lobe for both stars, filled grey if the star fills its Roche 
lobe. The ruler at the bottom shows the scale. 



We performed spectral fitting of both individual com- 
ponents and derived T c s — 10800 ± 550 K for the more 
massive component B and T c s = 8900 ± 550 K for compo- 
nent A. A ges ca n be estimated from the WD cooling tracks 
of Wood ( |1995|) . The resulting cooling ages for both white 
dwarfs are close to 1 Gyr. 

We compared the resulting orbital parameters of 
HE 1414-0848 with the outcome of theoretical cal- 
culations for the formation of DDs and find that 
HE 1414—0848 seems to be a fairly typical double white 
dwarf in terms the masses of the two white dwarfs and its 
orbital period. 
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